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ABSTRACT 

The primary goal of the Taiwan ECDFS Near-Infrared Survey (TENIS) is to find weU screened 
galaxy candidates at 2; > 7 {z' dropout) in the Extended Chandra Deep Field-South (ECDFS). To 
this end, TENIS provides relatively deep J and Kg data (~ 25.3 ABmag, Scr) for an area of 0.5 x 0.5 
degree. Leveraged with existing data at mid-infrared to optical wavelengths, this allows us to screen 
for the most luminous high-z objects, which are rare and thus require a survey over a large field to be 
found. We introduce new color selection criteria to select a z > 7 sample with minimal contaminations 
from low-z galaxies and Galactic cool stars; to reduce confusion in the relatively low angular resolution 
IRAC images, we introduce a novel deconvolution method to measure the IRAC fluxes of individual 
sources. Illustrating perhaps the effectiveness at which we screen out interlopers, we find only one 
z > 7 candidate, TENIS-ZDl. The candidate has a weighted Zphot of 7.8, and its colors and luminosity 
indicate a young (45M years old) starburst galaxy with a stellar mass of 3.2 x 10^° M©. The result 
matches with the observational luminosity function analysis and the semi-analytic simulation result 
based on the Millennium Simulations, which may over predict the volume density for high-z massive 
galaxies. The existence of TENIS-ZDl, if confirmed spectroscopically to be at z > 7, therefore poses 
a challenge to current theoretical models for how so much mass can accumulate in a galaxy at such a 
high redshift. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Finding objects at high redshifts is of fundamental im- 
portance in observational cosmology. By doing so, we 
are probing closer to the first galaxies and black holes in 
the universe, and therefore ever closer to the initial condi- 
tions for the formation of galaxies and structures. At the 
present time, we do not know when galaxies first formed, 
nor how quickly they can grow in mass during their ex- 
treme youth. Because z 7 is the epoch that marks 
the end of the reionization, finding objects at z > 7 will 
provide crucial constraints on the source of ionizing pho- 
tons, as well as in the early assembly history of galaxies 
and black holes. 

The last decade or so has bore witness to an explosive 
growth in the discovery of galaxies or galaxy candidates 
at z ^ 6 — 10, pushing our understanding of galaxy 
formation to ever earlier epochs after the Big Bang. 
Many efforts have been made using various techniques 
(e.g., iTaniguchi et al. I l2005t iKashikawa et al. I 120061 : 

^ Institute of Astrophysics & Astronomy, Academia Sinica, 
P.O. Box 23-141, Taipei 106, Taiwan, R.O.C. 

^ Department of Physics and Astronomy, University of Mis- 
souri, Columbia, MO 65211, USA 

^ National Astronomical Observatory of Japan, 2-21-1 Osawa, 
Mitaka, Tokyo, 181-8588, Japan 

^ Institute for the Physics and Mathematics of the Uni- 
verse, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa,277- 
8583 Japan 

^ Department of Physics, University of Hong Kong, Pokfulam 
Road, Hong Kong 

® Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138, USA 

Infrared Processing and Analysis Center, Cahfornia Institute 
of Technology. 770 South Wilson Avenue. Pasadena, CA 91125, 
USA 



Ive et al. I [20061: iShimasaku et al. I 120061: 
2007: 'Mannucci et al. I '2007': 'Ota et al. I 



Richard et a L,, 20081: Bouwcn s et a' 



2010s 



Stark et al. I 
200l 



al lbllc]: lOesch et al I [ 2009. 201013 



20101: iMcLure et al I 12009., ,2010n20 
MalHickev et al. 11201 Ot iHenrv et al. 



[2oa 



2010; 



2009 



Castellano et al. 
^ Bunker et al. 
2007L 12008. 200^ 



Bradlev et al. ' '2008': 'Zheng et al. "20091: iCapak et al 



2011; Wilkins et al. 2010, 2011; Sob ral et al. T 1 2009 : 
Gonzalez et al. I 120101: lOuchi et al. I 120091 I2OIO : 
Yan et al. I 1201 If ), but all can be roughly separated 
into three categories: (1) ground-based narrow-band 
imaging searches for Lya emitters (LAEs) (2) space- 
based (HST/WFC3) ultra-deep narrow-field surveys 
for Lyman Break galaxies (LBGs) (3) ground-based 
deep and wide-field surveys for LBGs. The continuum 
emissions of the LAEs provided by category (1) are 
usually too faint to be detected even photometrically; 
and so their UV star-formation rate and the continuum 
luminosity function (LF) of these objects cannot be 
derived. After HST/WFC3 was launched, the sample 
size of category (2) increased rapidly and now dominates 
the z > 6 sample. Given the limited survey field 
sizes, however, these studies mainly focus on very faint 
samples as the probability of finding luminous but rare 
objects are small. By contrast, surveys corresponding to 
category (3) are capable of finding rare bright sources 
that provide the strongest leverage in constraining the 
bright end of the LF. 

Most of the current z > 7 sample are provided by sur- 
veys in category (2) and (3). The vast majority of these 
z > 7 candidates have been identified based on their 
photometric colors and do not have spectroscopic con- 
firmations. Contamination by intervening objects may 
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be severe in the currently known sample. There are two 
major sources of contamination, heavily reddened low-^ 
galaxies and Galactic cool stars. For the samples pro- 
vided by category (2), low-z galaxies can be well resolved 
by HST/WFC3 and so their identification is straightfor- 
ward. A Galactic cool star, however, may have both 
similar colors and morphology as a compact galaxy at 
z > 7, and so contamination from cool stars remains an 
issue. Samples provided by category (3) suffer from these 
major contaminations. This issue ha s been d iscussed in 
detail in many prev ious studies (e.g., |Ouclii_et al. 2009; 
ICapak et al. 1120 111) , and has been recognized to be diffi- 
cult to resolve. In this paper, we use Spitzer IRAC pho- 
tometry to introduce additional criteria for significantly 
reducing both the abovementioned contaminations in the 
category (3) studies for selecting z > 7 candidates. 

As we utilize the Spitzer IRAC photometry in our 
analysis, we require a field with deep multi-wavelength 
data including Spitzer IRAC observations. For this pur- 
pose, the Exte nded Chandra Deep Field South (ECDFS, 
iLehmer et aLl bOOS) is ideal. The ECDFS is a 0.5 x 0.5 
degree field that surrounds both the Chandra Deep Field 
South (CDFS) and the Great Observatories Origins Deep 
Survey South (GOODS-S) field. The ECDFS has been 
studied at many wavelengths in addition to the origi- 
nal Chandra X-ray data. In the UV, the Deep Imaging 
Survey (DIS) for the Gal axy Evolution Explorer ultra- 
violet sateUite (GALEX, iMartin et ahl [2005) observed 
the entire ECDFS. In the optical, the Galaxy Evolution 
from Morphology an d SEDs (GEMS, .Rix et al. . ,200l 
iCaldwell et al. 1120081) project used the Advanced Cam- 
era for Surveys (ACS) on the Hubble Space Telescope 
(HST) to cover nearly the entire ECDFS. The ECDFS 
has also been obserbed in many ground-based optical and 
near-infrared imaging projects, e.g., the Classifying Ob- 
jects through Medium-Band Observati ons, a spectropho- 
tometric 17-fiher survey fCOMBO-17. lWolf et aLl[200ll 
optical only) pr ogram, Multiwayelength Survey by Yale- 
Chile (MUSYC.IGawiser et al. l[2006l:TTavlor et al. l[200l 
iCardamone et al. I l2010( ) including 10 0/NIR broad- 
bands and 18 deep topical medium-bands, and the NIR 
survey using ISAAC on VLT in the GOODS-S field. 
In the infrared and far-infrared, most of the relevant 
data have been obtained using Spitzer. The Spitzer 
IRAC/M USYC Public Legac y in the ECDFS survey 
(SIMPLE. Da men et al.1 l2011|) provides deep IRAC data 
at wavelengths of 3.6, 4.5, 5.8, and S.Ofim with a 
10' X 15' ultra-deep GOOD S-S IRAC data in the cen- 
ter [Dickinson et al. I ()2003D , and the F ar-Infrared Deep 
Extragalactic Legacy survey (FIDEL, [Dickinson et al. I 
I2007D took MIPS images at wavelengths of 24, 70, and 
160/zm to cover nearly the entire ECDFS. At radio 
wavelengths, th e ECDFS has been cove red by the VLA 
1.4GHz Survey (jKellerman et al. l[200l . 

Like in previous surveys in category (3), we start by 
using the Lyman-break technique to find LBGs at z > 7. 
The wavelength of the Lyman-break is shifted to between 
z band and J band for LBGs at z = 7 — 10, i.e., LBGs 
at z = 7 — 10 have very red z — J colors (z — J > 2, 
all the magnitudes and colors in this paper are in AB 
system, unless noted otherwise). Hence z-dropout ob- 
jects are what we look for. In the ECDFS, the deep- 
est available z-band data is the HST/ACS_F850LPdata 
from the GEMS (,Rix et al. ..2004: .Caldwell et al. |[2008l) . 



which has a 5a limiting magnitude of 27.1. To gain the 
full advantage of the deep F850LP data, J-band data 
with a 5fT limiting magnitude of 25.1 are a minimal re- 
quirement. In the ECDFS, two J-band data samples are 
available: the GOODS-S data taken using ISAAC on the 
VLT, and the MUSYC data. The depth of the GOODS-S 
J band is about 25.1 mag, which meets our requirement. 
The field size of the GOODS-S, however, is only about 
one- fifth of that of the ECDFS, significantly reducing the 
probability of finding rare and luminous LBGs by a fac- 
tor of five. On the other hand, although the J-band data 
provided by the MUSYC covers the entire ECDFS, the 
depth is only about 23.0 mag, and thus far too shallow 
for finding z > 7 objects. A new deep (J > 25.1) near- 
infrared imaging survey with the full ECDFS coverage is 
therefore needed for our study. 

In this paper, we describe the deep near-infrared survey 
that we carried out in the ECDFS, and the data reduc- 
tions performed on the optical, near-infrared, and IRAC 
data, in § [2] In § [Sj we provide a description of the sam- 
pling criteria for this study. The possible contaminations 
in our sample are discussed in § |4l We present our re- 
sults and the implications of our sample in § [5l In § |6l we 
summarize our results and discuss future work. The cos- 
mological parameters used in this study are fl^ = 0.73, 
Qm = 0.27, Ho = 71 km s^^ Mpc-\ and w = -1. 

2. DATA AND PHOTOMETRY 

2.1. TENIS Data 

We initiated the Taiwan ECDFS Near-Infrared Survey 
(TENIS) in 2007. The goal of this project is to obtain 
deep J and Kg images in the ECDFS, and thus fill in the 
big wavelength gap between 1.0/xm and 3.6/im. TENIS 
is yet the deepest near-IR survey in a 30' x 30' sky area 
with a full coverage of all four IRAC channels. Figure [T] 
shows the fields of the main surveys in the ECDFS: SIM- 
PLE, GOODS-S(ACS and IRAC), MUSYC, COMBO- 
17, GEMS, the TENIS J, K„ and Y. The scientific goal 
of the TENIS J-band observation is to search for the 
most luminous high-z objects (i.e., this work), which are 
rare and thus require a survey over a large field to be 
found, while the TENIS if^-band observation is for the 
studies of the mass assembly history and sub-millimeter 
galaxies at 2 > 2. After combined with existing deep 
images in the optical and 3.6-8/xm, the TENIS data can 
be used to derive robust photometric redshifts and stel- 
lar masses of a large sample of galaxies. Therefore the 
TENIS dataset is unique not only in its depth and survey 
area but also in its rich multi-wavelength ancillary data 
in the ECDFS. 

2.1.1. CFHT WIRCam Observations 

The TENIS data were take n using the Wide- field 
InfraRed Camera (WIRCam, iPuget et al.1 I200I on 
the Canada- France-Hawaii Telescope (CFHT). WIRCam 
consists of four 2048x2048 HAWAII2-RG detectors cov- 
ering field of view of 20' x 20' with a 0".3 pixel scale. 
To date we have spent 44 hours in J and 40 hours in 
Ks. The J-band data were taken in 2007B and 2008B, 
with an average 5cr limiting magnitude of J = 25.5 for 
point sources. The Xs-band data were taken in 2009B 
and 2010B, with an average 5(t limiting magnitude of 
Kg — 25.0. The average seeings were 0".8 and 0".7 for 
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Fig. 1.— The TENIS J-band image. The black lines indicate the TENIS J and X^-band field. The green lines delineate the TENIS 
y-band field. The blue, pink, cyan Unes indicate the fields of GEMS, COMBO-17, and MUSYC, respectively, and the GOODS-S (IRAC 
and ACS) fields are represented by brown and dark green lines. 

J and Kg, respectively. More details of the WIRCam 
imaging observations is described in the TENIS cata- 
logue paper (Hsieh et al., in preparation). 

2.1.2. Subaru MOIRCS Observation 

We obtained Y data in the ECDFS using the Subaru 
MOIRCS in 09B. The observation is a follow-up observa- 
tion of the TENIS project. The goal of the observation 
is to have F-band photometry for a set of high-z objects 
pre-selected using their z — J colors. Two pointings were 
planned to be observed. One is inside the GOODS-S 
field and the other is outside the GOODS-S field. How- 
ever, due to the bad weather, only one pointing on the 
east side of the ECDFS was observed, i.e., a 4' x 7' area 
outside the GOODS-S field. The 5a limiting magnitude 
for point sources is 25.5. Although the Y data are only 



available for part of the ECDFS, the data are still very 
important to further improve the photometric redshift 
quality for the candidates in that area since the wave- 
length of Y band is right between those of F850LP and 
J. 

2.1.3. Data Reduction 

The TENIS J and Kg data were processed using an In- 
teractive Data Language (IDL) based reduction pipeline 
called Simple Imaging and Mosaicki ng Pipeline. T he de- 
tails of the pipeline are described in ' Wang et al. I ([2010). 
We used the pipeline to deal with fiat-fielding, removing 
instrumental features like crosstalk and residual images 
from saturated objects in previous exposures, masking 
satellite trails, distortion and astrometry correction, mo- 
saicking and stacking images, and photometric calibra- 
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tion. The zero-point calibration was done by matching 
the fluxes in an aperture of 5" in diameter with those in 
the Two Micron All Sky Survey (2MASS) point source 
catalogs. 

2.1.4. Photometry 

We used SExtractor version 2.5.0 ()Bertin fc Arnoutsi 
119961 to detect objects and measure their fluxes in the 
WIRCam images. For the TENIS J image, we used the 
"FLUX_AUTO" values of the SExtractor output as the 
measured J fluxes. The flux errors provided by SEx- 
tractor do not include the correlated errors. Therefore 
we calibrated the errors with the following procedure: 
First, the fluxes and the flux errors were re-measured us- 
ing SExtractor with a 2" diameter aperture. We then 
convolved the source-masked image with a 2" diame- 
ter aperture, and calculated the rms around each pixel 
on the convolved images. The ratio between the aper- 
ture photometric flux error provided by SExtractor for 
a certain object and the rms value around the same po- 
sition in the convolved image is the correction factor of 
the flux error for that object. The median value of the 
correction factors was computed to be the general cor- 
rection factor for all the sources. For J-band data, the 
general correction factor is 0.793. The same procedures 
were done with different aperture sizes. We find that the 
correction factors are very stable with different aperture 
sizes, which is consistent with the experience in Wang 
et al. 2010. Therefore we just applied the factor to the 
"FLUXERR.AUTO" values for all the TENIS J objects 
to calibrate their flux errors. 

For the TENIS Kg image, the double-image mode 
of SExtractor was performed to generate a J-selected 
Ks photometric catalogue. In order to use the double- 
image mode, we shifted the Kg image to match to the 
J image. The Ks fluxes were also measured with the 
"FLUX_AUTO" option, and the flux errors were cali- 
brated with the same procedure used in the TENIS J 
data. 

To estimate the completeness for our data, we used a 
Monte-Carlo simulation. The faint objects are compact 
but may not be all point-like. The medium full width 
half maximum (FWHMs) of objects with J 25 is about 
2".0 (~ 10 kpc at z ~ 7). Therefore, to better determine 
the completeness of faint objects, we randomly put 100 
artificial sources into the TENIS J image, with an im- 
age that is averaged from ~ 100 real extended sources 
with FWHMs between 2".0 and 2". 2, and with signal- 
to-noise ratios greater than 100. The magnitudes of the 
artificial sources are uniformly distributed in the mag- 
nitude space from 20mag to 27mag, and Poisson errors 
were added into the images of the artificial sources. We 
then ran SExtractor to do object finding. This proce- 
dure was run for 1000 times to deliver the detection rate 
of these artificial sources in each magnitude bin, i.e., the 
completeness. Figure [2] shows that the surface density 
vs. apparent magnitude of the TENIS J data. The com- 
pleteness corrected surface density is used in our analyses 
later. 

2.2. Ancillary Data 

We incorporated data from several ECDFS surveys 
at different wavelengths to perform our analyses. To 
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Fig. 2. — The surface density vs. TENIS J-band magnitude plot. 
The surface density is measured by selecting sources using SEx- 
tractor. The large open box indicates the surface density for each 
magnitude bin. The small open box indicates the surface density 
for each magnitude bin after the completeness correction. The dot- 
ted line is the best fit of the surface density after the completeness 
correction. 

measure accurate colors for the color selection, we pro- 
cessed the images from the ancillary data to generate 
new images that are pixel-to-pixel matched to the TE- 
NIS images. For the space-based optical images, we also 
matched their PSFs to the TENIS data to ensure better 
color measurements. In this subsection, we describe how 
we processed the ancillary data and measured fluxes. 

2.2.1. MUSYC Data 

The MUSYC team has released their deep optical and 
near- infrared images in the ECDFS. To improve the qual- 
ity of the photometric redshifts, which will help eliminate 
foreground contaminations in our study, we included the 
MUSYC optical data in our analyses. We resampled the 
MUSYC images in f/, f738, B, V, i?, /, z, and O3 to 
match the TENIS J image, and then used the double- 
image mode of SExtractor with "FLUX_AUTO" to gen- 
erate a TENIS J-selected MUSYC photometric cata- 
logue. The zero-point calibration was done by match- 
ing the fluxes with t hose in the pubhshed BVR-selecied 
MUSYC catalogue (jCawiser et al. II200I ). The flux er- 
rors of the MUSYC photometry were calibrated using 
the method described in § 12.1.41 

2.2.2. GEMS ACS Data 

Th e GEMS project (|Rix et al. l[200l ICaldweU et aD 
I2008f) provides deep HST ACS images with the F606W 
and F850LP fllters in the ECDFS, except some small 
areas close to bright objects. The 5a limiting magni- 
tudes for point sources are 28.4 and 27.1 for F606W^ and 
F850LP, respectively. We mosaicked the public GEMS 
ACS images and then resampled the images to match the 
TENIS J image. The PSFs of the ACS images are very 
small compared to the PSFs of the WIRCam images. To 
derive better color measurements between the ACS fll- 
ters and J, we applied a smoothing kernel to the ACS 
images to match their PSFs to that of the TENIS J im- 
age. We then repeated the same procedure as what we 
did for the MUSYC data to deliver the J-selected GEMS 
photometric catalogue with calibrated flux errors. 
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2.2.3. SIMPLE IRAC Data 

The SIMPLE project (jPamen et al. 1 120 111 ) provides 
deep IRAC observations covering the entire ECDFS with 
the 10' X 15' G OODS-S IRAC mosaics in its center 
iDickinson et al. I p003). We mosaicked the SIMPLE 
images and resampled them to match the TENIS J 
image. For the flux measurement, however, we could 
not use SExtractor in the double-image mode with the 
"FLUX_AUTO" option to generate the J-selected IRAC 
photometric catalogue, since the PSFs of IRAC images 
are about 1".5 to 2".0 arcsec, which are much larger than 
that of the TENIS J image. If we were to do this, the 
IRAC fluxes would be seriously under-estimated as SEx- 
tractor uses the aperture size estimated from the refer- 
ence image (TENIS J image) to measure the fluxes on 
the target image (SIMPLE IRAC images). 

One can use SExtractor in the double-image mode with 
the standard aperture photometry to measure the IRAC 
fluxes, and then apply the aperture correction values pro- 
vided in the IRAC data handbook (Chapter 5, Table 5.7) 
to derive the estimated IRAC total fluxes. However, due 
to the very large PSFs of IRAC images, using simple 
aperture photometry to measure the fluxes of deep IRAC 
images in a relatively crowded field like the ECDFS has 
serious confusion issues and the local background estima- 
tion used in aperture photometry can be easily affected 
by neighboring objects. We note that the accuracy of the 
IRAC photometry is critical in our analyses (see § [3|), 
requiring an alternative method to estimate the IRAC 
fluxes properly. 

Many intensive efforts have gone into developing meth- 
ods for estimating fluxes in an image with a large PSF 
(imagel, hereafter). The most recent method relies on 
utilizing the morphological information from a high qual- 
ity image in anothe r band (imaKe2, hereafter), (e.g., 
Grazian et aTl l2006t iLaidler et al. 1 l2007t iWang et all 
20101 : IMcLure et al. 1 I2011D . Bv assuming the intrinsic 
morphology of objects are identical in the two wave- 
lengths (i.e., no color gradient from its center to edge), 
one is able to use the PSF of image2 to deconvolve that 
image, and then convolve the deconvolved image2 with 
the PSF of imagel. The scaling factor that needs to be 
applied to match the peak values between an object in 
image2 and the same object in the processed imagel is 
the flux ratio (color) of that object. The assumption 
that there is no color gradient, however, does not always 
stand. We have therefore developed a new method that 
is quite different from what has been used, but is very 
similar to the traditional CLEAN deconvolution in radio 
imaging. 

We used the TENIS J image as a prior image (i.e., 
image2) to estimate the fluxes in the IRAC images (i.e., 
imagel). The segmentation map of the TENIS J image 
generated by SExtractor with a 0.7(7 threshold preserves 
the information of what pixels are attributed to each ob- 
ject, and it provides the positional and shape/boundary 
information of each object when deconvolving the IRAC 
images. We only used the boundary information rather 
than the full morphological information adopted in pre- 
vious works. The boundary information introduced with 
the segmentation map is often called the "CLEAN win- 
dow" in radio image. The IRAC PSFs were generated 
using bright isolated IRAC point sources in each IRAC 



channel. The size of the PSF images is 1' x 1'. Due 
to the non-spherical IRAC PSF and different orienta- 
tions of Spitzer in the observing epochs for GOODS-S 
and SIMPLE, the PSFs of the stacked IRAC image in 
GOODS-S region is very different from those in the ex- 
tended area. We tried to use the PSF generated from 
the SIMPLE data to deconvolve the GOODS-S region 
and found that the fluctuation of the residual image (see 
the next paragraph for the detail of the residual image) 
is 5 times larger (i.e., ~ 1.7 mag worse in limiting magni- 
tude) as compared to the result using the PSF generated 
from the GOODS-S data itself. We therefore generated 
an individual PSF for each region and deconvolved the 
two regions separatedly. 

The deconvolution process always starts at the pixel 
with the highest value measured within a 9 pixel aper- 
ture in diameter {Fap9^ hereafter), and this pixel must 
have been registered to one object in the J segmenta- 
tion map. Since the SIMPLE image (~ 0".6 per pixel) 
has been resampled to match the TENIS J-band image 
(~ 0".3 per pixel), simply moving a window (i.e., FApg) 
across the SIMPLE image to find the location of a peak is 
very similar to do a sub-pixel centering. Once this pixel is 
found, we subtracted a scaled PSF from the surrounding 
1' X 1' area centered on this pixel. The ratio of the sub- 
tracted fiux depends on Fap9] if Fap9 is greater than 5a, 
then 0.1% of the flux was subtracted. If Fap9 is less than 
5(7, then 100% of the flux was subtracted. The subtrac- 
tion ratio that we use (0.1%) is much smaller than the 
values that most deconvolution/clean algorithms usually 
use (10%). We found that a subtraction ratio of 10% 
for our deconvolution process sometimes produces much 
worse results especially for bright extended sources be- 
cause our algorithm is less restricted; 0.1% is the value 
that has a good balance between performance and qual- 
ity. The subtracted flux (Fsub) was summed and reg- 
istered to a certain reference object according to the J 
segmentation map. After the subtraction, the deconvo- 
lution process is repeated on the subtracted image, until 
there was no pixel with Fap9 higher than 1.5ct. We note 
that the concept of the deconvolution here is identical 
to that of CLEAN in radio imaging. At the end of the 
deconvolution process, the Fsub for each object is the 
flux measurement of that object, and the flnal subtracted 
image is the residual map. The latter allows us to check 
the quality of the deconvolution, and it is also used for 
estimating accurate photometric errors. The flux error 
of each object was calcualted based on the fluctuation of 
the local area around that object in the residual map. 
Any imperfection of the PSF would cause larger fluctua- 
tion in the residual map, and that was taken into account 
by the flux error calculation. 

We demonstrate the performance of this method in 
Figure |3l A small region of the TENIS J image, the 
SIMPLE 3.6^m image, and the residual SIMPLE 3.6^m 
image are shown in the left, middle, and right panels. 
The brightness and contrast of the three panels are ex- 
actly the same. By comparing the SIMPLE 3.6/im image 
to the TENIS J image, one can see that many objects 
in the most dense area are mixed together in the 3.6/im 
image while they are still well-separated in the TENIS 
J image. After processing the 3.6/im image using the 
method described above, the residual image as shown in 
the right panel of Figure [3] shows that the deconvolution 
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works very well. Sources left in the residual image are 
objects detected in the IRAC images but not the TENIS 
J image. An example is the object sitting at the center 
of the red circle in the right panel of Figure |3l This is a 
bonus of the deconvolution method. 

Figure |4] shows a comparison of IRAC photometry us- 
ing the abovementioned deconvolution method against 
that using the SExtractor with FLUX_AUTO. As can be 
seen, the results of two photometric methods agree with 
each other for bright objects, but the scatter between the 
two measurements increases for fainter objects. This is 
because the flux measurements of faint objects are much 
more easily affected by neighboring objects, coupled with 
the fact that fainter objects have higher mean surface 
number densities. Some bright sources have > 0.1 mag 
differences between their deconvolution and auto magni- 
tudes. According to the flags provided by SExtractor, 
more than 80% of objects brighter than 20 mag are se- 
riously blended with their neighbors. We also confirmed 
that all the objects having large differences between their 
deconvolution and auto magnitudes are blended with 
their neighbors as shown in their flags. The deconvo- 
lution method can deliver more accurate flux estimation 
for faint sources in the IRAC images, an important fac- 
tor in our selection criterion for z > 7 objects (see §|31)- 
A more detailed description of the deconvolution method 
we used to estimate the IRAC fluxes is described in Hsieh 
et al. (in preparation). 

2.3. J-selected master catalogue 

We combined all the data described in this section and 
generated a 16-band J-selected master catalogue con- 
taining 114K objects. The 5cr limiting magnitude for 
point sources in each band is shown in Table [TJ For the 
MUSYC, GEMS, and TENIS data, the values were de- 
rived using an 1".0 aperture in diameter wit h aperture 
and n oise correlati on correc t ion. A ccording to lRix et al. I 
( |2004D and CaldwelTeFaD ()2008[ ). the 5cr limiting mag- 
nitudes for the GEMS F606VK and F850LP are 28.53 
and Tl .21 ^ which are much deeper than those listed in Ta- 
ble [1] because the aperture size that they use is ^ 0".l 
in diameter, which is 10 times smaller as compared to 
ours. For the SIMPLE data, the aperture we used in the 
IRAC deconvolution process is a 9-pixel window, which 
mimics an 3".0 aperture in diameter for point sources. 
The IRAC limiting magnitudes in Table [1] which are the 
medium values derived using point sources with 5(t detec- 
tions in both SIMPLE and GOODS-S IR AC data. The 
values are different from those listed in iDamen et al. I 
(|201lD because of the different methodologies. More de- 
tails are provided in the TENIS catalogue paper (Hsieh 
et al., in preparation). This master catalogue is used to 
select our z > 7 sample. 

3. CANDIDATE SELECTION 

3.1. Color criterion 

For z-dropout studies, the most important issue is to 
eliminate contaminations from low-z galaxies and Galac- 
tic cool stars. Most previous studies make great efforts to 
estimate the contamination rates, which are not negligi- 
ble. The accuracies and uncertainties of these estimates, 
however, are unknown; the best way to deal with the 
contaminations is to eliminate them in the first place as 



TABLE 1 

5(T LIMITING MAGNITUDES FOR POINT SOURCES 
IN THE J-SELECTED MASTER CATALOGUE 



Band 5cr limiting magnitude 

MUSYC U 26.44'' 

MUSYC U38 26.34'' 

MUSYC B 27.21'' 

MUSYC V 26.85" 

MUSYC R 26.80" 

MUSYC I 25.06" 

MUSYC z 24.29" 

MUSYC OS 25.88" 

GEMS ACS F606W 26.80*' 

GEMS ACS F850LP 25. 72^ 

TENIS J 25.57 

TENIS Ks 24.95 

IRAC 3.6iim 25.82 

IRAC 4.5fj.m 25.75 

IRAC 5.8fj,m 23.77 

IRAC S.Oirni 23.62 



" The MUSYC limiting magnit udes in this ta- 
ble ar e de eper than those in IGawiser et al. I 
SWOB) and lTavlor et al. 1 1|2009|) because we use 
a smaller aperture size. 

^ The GEMS limiting magni tudes in th i s table 
are shallower than t hose i n IRix et al. I I I2004I) 
and ICaldwell et al. I l|2008l ) because we use a 
much larger aperture size. See text for details. 

far as is possible. The IRAC photometry provide a good 
solution to this issue and we discuss how to utilize the 
IRAC photometry to select a relatively pure sample in 
this subsection. 

Galaxies at different redshifts in a J — 3.6/im vs. 
F850LP — J diagra m are shown in Figure [5l The 
GALAXEV template fBruz ual fc Chariot N2003I ) is used 
to generate this plot. LBGs at z > 5, late-type, and 
early-type galaxies at z < 5 with different extinctions 
are plotted. Most of the previous z-dropout studies use 
a simple color cut (e.g., z — J > 2.0) to select z > 7 ob- 
jects. However, as shown in Figure [Sj heavily reddened 
low-z galaxies would have very red z — J colors, and the 
z — J color of a late-type galaxy with Av > 3 would 
be very similar to that of z > 7 objects {z — J > 2). 
Our two-color criteria in the J — 3.6/im and F85QLP— J 
space, however, greatly minimizes this issue. The selec- 
tion criterion are accordingly 

F850LP - J > 2.0, 
J — 3.6/iTO > —1.0, 
F850LP- J > 0.8 X (J- 3.6/im) + 1.1, (1) 

as shown by the dark gray area in Figure [5] The third 
criterion is closely parallel to the reddening direction of 
low-z galaxies, and therefore greatly reduces the contam- 
ination rate from low-z galaxies. 

Although contamination by low-z galaxies can be ef- 
fectly dealt with using the criteria defined in Equation [1] 
contamination by Galactic cool stars can still be prob- 
lematic. We use the AMES-dusty dwarf star template 
{AUard et al. 2001) to estimate the distribution of M, L, 
and T dwarfs in Figure [5l which is plotted as a light gray 
area. The observed M and L dwarfs from IRTF Library 
are also shown as dark star signs. We note that there is 

^ |http : //irtf web . if a . hawall ■ edu/$\slm$spex/IRTF_Spectral_Llbrary71 
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Fig. 3. — These three panels show the deconvolution quaUty of the IRAC images. A small region of the TENIS J image, the SIMPLE 
image, and the residual SIMPLE image are shown in the left, middle, and right panels. The region of each panel is 135" X 135". 

North is up and east is to the left. The brightness and contrast of the three panels are exactly the same. The clean residual image in 
the right panel shows that the deconvolution works very well. The unsubtracted object marked by the red circle in the right panel is an 
IRAC-bright source but that is not detected in J. 
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Fig. 4. — The comparison of the IRAC photometry using differ- 
ence methods. The differences between the magnitudes using the 
deconvolution method (decon) and that using the SExtractor with 
FLUX_AUTO (auto) are shown. It shows that the photometry us- 
ing two different methods are consistent for bright objects but the 
scatter between the two measurements increases for fainter objects. 



an overlapped region between our selection criteria (dark 
gray area) and the dwarf distribution (light gray area). 
According to the AMES-dusty template, dwarfs in the 
overlapped region are L dwarfs with effective tempera- 
tures between 1700K and 1900K. Therefore, we need to 
make another criterion to discriminate against very cool 
dwarf stars. 

Fortunately, the IRAC colors of dwarf stars are very 
different from those of high-z galaxies. Therefore we 
can minimize the contamination issue from Galactic cool 
stars by adding another selection criterion: 

3.6/iTO - 5.8nm > 0.12 x ( J - 3.6fim) - 0.4. (2) 

Figure ini shows how this additional criterion works. The 
solid and dashed lines indicate starburst galaxies at z — 
5 — 10 with Ay = and Ay = 1, respectively. The 
arrow indicates the reddening direction. The dark gray 
area shows the selection criteria for high-z galaxies while 
the light gray area shows the region occu pied by dwarf 
stars as computed using the AMES-dusty (|Allard et al. I 
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Fig. 5. — Selection criterion in J — S.dfim vs. F850LP — J dia- 
gram. The thick solid line indicates LBGs at 2 > 5.0 and the thick 
dotted line indicates LBGs with Av = 1 at z > 5.0 The two thin 
solid lines indicates early-type and late-type galaxies a.t z < 5.0, 
and the thin dotted line indicates late-type galaxies with Av = 1 at 
z < 5.0. The black arrow shows the reddening direction. The dark 
gray area represents the selection criteria of Equation [T] while the 
light gray area indicates where the Galactic cool stars are, whic h 
are estimated using the AMES-dusty model HAUard et al. ||2001| V 
The observed M and L dwarfs from IRTF Library are plotted with 
dark star signs. The red dot indicates the candidate that we find 
(see §|4]for the details). This figure shows that although the selec- 
tion criteria can pick up LBGs at z > 7 and avoid low-z galaxies 
with any level of reddening, some L dwarf stars could still be se- 
lected as well. 



120011 ) model. In FigurelHl high-z galaxies and dwarf stars 
are well separated, demonstrating that IRAC colors can 
break the color degeneracy between high-z galaxies and 
dwarf stars in the optical and near-infrared bands. 

To improve the quality of the candidate selection, addi- 
tional criteria are added to eliminate low-z galaxies that 
satisfy the above selection criterion because of photomet- 
ric errors. First, all candidates must have more than Su 
detection in the WIRCam J band. Second, all candi- 
dates must have fluxes lower than l.Scr in optical bands 
bluer than F850LP. We also stack all the optical images 
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Fig. 6. — IRAC color-color diagram for high-z galaxies and dwarf 
stars. The solid line and dashed line indicate starburst galaxies at 
z = 5 — 10 with Ay = and Ay = 1, respectively. The arrow 
indicates the reddening direction. The dark gray area shows the 
selection criteria for high-z galaxies while the light gray area shows 
where the dwarf stars are. The red dot indicates the candidate that 
we find (see § 3] for the details). According to this figure, using 
IRAC colors can successfully break the color degeneracy between 
high-2 galaxies and dwarfs in optical and near infrared bands. 

(except the F850LP and MUSYC z bands as objects 
at z > 7 may still be detect ed in z bands) to gen erate 
an ultra-deep optical image. 'Mc Lure et all (|2011[ ) have 
found that some previously published z-dropout candi- 
dates actually have counterparts in ultra-deep stacked 
optical images, that are just contaminations. The ultra- 
deep stacked optical image therfore allows us to check 
if some contaminations cannot be detected in our sin- 
gle optical bands because of insufficient sensitivity. We 
generate two different stacked optical images: one is av- 
eraged from MUSYC C/, t738, B, V, R, I, and GEMS ACS 
FQQQW (7 filters); the other is averaged from just red 
filters (MUSYC and GEMS ACS F606VF). The 

former provides the deepest stacked optical image while 
the latter is optimized for detecting faint red objects. At 
this stage, 143 objects satisfy all the color criteria. 

3.2. Photometric redshift criterion 

Although the two color-color selection criteria are very 
robust for picking up z > 7 objects, contaminations can 
still be present because of photometric errors and tran- 
sients. Photometric redshift can help to further reduce 
these issues. It utilizes all the data (16 bands) rather 
than just the 4 filters used in our color-color selections, 
so it is less affected by ph otometric errors and tr ansients. 
We use the EAZY code (jBrammer et al. Il2008i^ to esti- 
mate the photometric redshift of our candidates. We 
allow a 5% systematic flux error for each band to take 
into account the systematic errors caused by the fact 
that for a given object we are not missing the same frac- 
tion of light in different passbands that have very dif- 
ferent spatial resolutions. To minimize the bias effect of 
the photometric redshift estimation due to different fit- 
ting templates, we use three different templates to run 
the EAZY code, thus obtaining three photometric red- 
shifts for each candidate. The three templates used are 



CWW-hKIN, defauh EAZY vl.O, and PEGASE. The 
CWW+KIN template is the CWW empirical template 
(jColeman. Wu. &: Weedmani 119801) with the extension 
prescribed bv iKinnev et al. I (|1996[ ). The default EAZY 
vl.O template is generated using the iBlanton fc Roweis I 
()2007[) algorithm with the PEGASE models and cali- 
brated using semi-analytic models, plus an additional 
youn g and dusty template. The PEG ASE template is 
fr omlFioc fc Rocca-Volmerange I ()1997D . 

lYan et al. I (|2004[ ) have found that many IRAC- 
selected extremely red objects (lERO) with very high 
photometric redshifts are actually galaxies at z < 2 pos- 
sessing two different stellar populations. The EAZY code 
can fit data with multiple stellar populations simultane- 
ously. For an lERO or similar object, using this function 
helps to break the degeneracy of high and low redshifts 
in the template fitting. We use this multiple popula- 
tion fitting function when employing the CWW-I-KIN 
and default EAZY vl.O templates. The documentation 
on the PEGASE template suggests that this template 
works better in the single population fitting mode, and 
so we decided not to use the multiple population fitting 
function with the PEGASE template. 

We compile a new list of candidates according to their 
photometric redshifts. We pick up an object with a pho- 
tometric redshift between 6.0 and 9.0, and the lower-limit 
of its photometric redshift at 68% confidence level must 
be greater than 3.0. Since there are three different pho- 
tometric redshifts for each object, all three photometric 
redshifts for a chosen candidate must meet the selection 
criterion simultaneously. We then assign the photomet- 
ric redshift with the smallest best-fitted chi-square as the 
redshift of each candidate. Applying the photometric 
redshift criterion left 22 objects in our candidate list. 

3.3. Cleaning the Sample 

All the previous selection criteria rely entirely on 
photometry (including the photometric redshift estima- 
tion). Photometric measurements, however, can be eas- 
ily affected by CCD bleeding, satellite tracks, diffraction 
spikes of bright stars, etc., which will also affect the per- 
formance of the selection criteria we use. To eliminate 
these false-detection cases, we visually check the images 
of every bands for all the 22 remaining candidates and re- 
move 2 affected objects from the list. In addition, objects 
detected in 24/im, X-ray, and 1.4GIIz are also omitted 
since they are very unlikely to be at high z. There are 3 
objects rejected at this stage. 

Although we use a color selection criterion in the J — 
3.6/im vs. 3.6fim — 5.8fim space to discriminate against 
Galactic cool stars, this criterion does not work well for 
some objects which are very faint in 3.6/irn and 5.8fim. 
This increases the contamination rate in the candidate 
selection because the uncertainties in their colors make 
them spill over into the selection region. To deal with this 
issue, we perform the SED fitting for the candidates using 
the AMES-dusty dwarf star template. We reason that 
an SED fitting would provide better information about 
the spectral type of an object as it uses the photometry 
from all the bands as compared to a pure color selection. 
We reject objects if they have similar or even smaller 
minimum chi squares with the AMES-dusty template as 
compared to those with the galaxy templates. 18 objects 
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are rejected in this way. 

In the end, only one object survives all the criteria im- 
posed. We name this object TENIS-ZDl. We show its 
thumbnails in Figure [7] and summarize its photometric 
results in Table [S] The listed detection limits for the 
optical bands are Icr values, and the size of the Kron 
aperture used for TENIS-ZDl is 3".2 x 1".7 (major and 
minor axes) according to the output of SExtractor. The 
FWHM of TENIS-ZDl in J provided by a Gaussian fit- 
ting is 2".3±0".46, demonstrating that it is an extended 
object (and hence not a star). The ultra-deep stacked 
optical images are also shown in Figure [5] to show that 
TENIS-ZDl is not detected in both ultra-deep optical 
images. As we mentioned in § 12.2.31 accurate IRAC 
flux measurements are essential for picking up our can- 
didates. Figure [7] shows that TENIS-ZDl has a bright 
close neighbor with a separation of just 3". Therefore the 
IRAC photometry of TENIS-ZDl measured using a typ- 
ical aperture photometry method would be seriously af- 
fected by this bright neighbor. In i) l2.2.3l we describe the 
deconvolution method that we use to derive IRAC fluxes, 
and Figure IH] shows a comparison between the original 
IRAC images around TENIS-ZDl, the IRAC images with 
TENIS-ZDl but without aU the J-band detected neigh- 
bors, and the IRAC residual images. It demonstrates 
that TENIS-ZDl is well-deconvolved in the IRAC im- 
ages and the photometric effect of its bright neighbor is 
minimized. The photometric redshift of TENIS-ZDl es- 
timated by EAZY is 7.822tJ:?^^, with a reduced x^=l-15. 
The best-fitted templates and the vs. redshift plots 
are sh own in Figure ITOl We also fit an AGN/quasar tem- 
plate CPol letta et al. II2007D and find a best-fit redshift of 
7.0691 

0.411- However, its is 6.50, which is higher than 
that of the galaxy fit. Although an AGN is not com- 
pletely ruled out, the fit favors a galaxy. 

4. CONTAMINATION 

In this sec tion we d i scuss the contaminations in our 
analysis. Ou chi et al. I (j2009t ) (Ouchi09, hereafter) lists 
all the possible sources of contamination for z-dropout 
studies. Here we follow the same items as in Ouchi09 to 
check the contaminations in our study. 

1) Spurious sources: One of the selection criterion of 
our z > 7 sample is that the objects must have more than 
3(7 detections in J, and TENIS-ZDl has a 4.3ct detection 
in J. Statistically, an object with a 4.3(7 detection is 
very unlikely a spurious source. As an additional check, 
we stack the images taken in 07B and 08B separately 
and generate two mosaic images for the two different 
semesters. We then measure fiuxes of objects in each 
image using the SExtractor in the double-image mode, 
by using the J-band master image (07B+08B) for the 
source detection. The magnitude of TENIS-ZDl in 07B 
is 25.03 ± 0.29, while that in 08B is 25.13 ± 0.24. Ac- 
cording to the fluxes of TENIS-ZDl in both semesters, 
it was well-detected in the two epochs. Furthermore, 
TENIS-ZDl has counterparts in Kg and IRAC images. 
All the above evidences suggest that TENIS-ZDl is not 
a spurious source. 

2) Transients: We repeat what we did for checking 
spurious sources. As we mentioned above, the flux of 
TENIS-ZDl in 07B is 25.03 ± 0.29, and that in 08B 
is 25.13 ± 0.24. The magnitude change between the 
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Fig. 7. — Thumbnail images for TENIS-ZDl. The region of each 
image is 10" X 10". North is up and east is to the left. TENIS- 
ZDl is not detected in all the optical bands and detected in all the 
WIRCam and IRAC images. 
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Fig. 8. — Ultra-deep stacked optical images for TENIS-ZDl. The 
left-hand image is averaged using MUSYC (/, f738, B, V, R, I, and 
GEMS ACS F606W images (7 filters). The right-hand image is 
combined from MUSYC V, R, I, and GEMS ACS F606VK images 
(4 red filters). The region of each image is 12" X 12". North is 
up and east is to the l eft. The details of the ultra-deep stacked 
optical images are in § 13.31 The red circle indicates the location 
of TENIS-ZDl in the J image. According to these two images, 
the optical flux of TENIS-ZDl is below the detection limit even in 
ultra-deep optical images. 

two epochs is only O.lmag, corresponding to a ^ 
10% flux variation. Furthermore, TENIS-ZDl did not 
move between the TENIS J-band (07B+08B), X^-band 
(09B-f lOB) and IRAC images, which shows that TENIS- 
ZDl is not a slow-moving solar system object. Therefore 
we conclude that TENIS-ZDl is not a transient. 

3) Low-z galaxies: Our selection criteria avoid select- 
ing galaxies at z < 6. However, colors of some low- 
z galaxies could also meet the selection criteria due to 
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TABLE 2 
Magnitudes of TENIS-ZDl 



U B R F606W F850LP J Ks 3.6tim A.5tim b.S/im S-O/im 

> 28.28 > 28.85 > 28.31 > 28.18 > 27.05 25. 12 ±.23 24.94 ±.34 24.29 ± .10 23.79 ± .06 23.28 ± .22 23.58 ± .40 
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Fig. 9. — The IRAC residual images of TENIS-ZDl. Images from left to right indicate images in IRAC 3.6fj,m, 4.5fim, 5.8fj,m, and 
S.Ofim. The upper images are the original images of IRAC channels. The middle ones are the images after subtracting all the neighbors 
(except the target) detected in J using the deconvolution method. The lower ones are the residual images after subtracting all the objects 
detected in J. The red circle indicates the location of TENIS-ZDl in the J image. The region of each image is 25" X 15". North is up and 
east is to the left. By comparing these three images, they show that TENIS-ZDl and its bright close neighbor are well-deconvolved and 
well-deblended, which suggests that the IRAC flux measurement effect of TENIS-ZDl due to the bright close neighbor is minimized. 

F850LP image. TENIS-ZDl is still far away from the 
area of low-z galaxies even if we take its error bars into 
account. Furthermore, although there is a local mini- 
mum at z < 
Figure [TOl its x 
than that at z 

type- 2 AGN at z ^ 1.5 can be very similar to that of a 
galaxy at z ~ 7, which makes it another possible source 
of low-z contaminants. We derived the expected fluxes at 
1.4GHz, 24^m, and 2-8keV for TENIS-ZDl by assuming 
a type-2 AGN at z 1.5; they are ~ 1785uJy (radio-loud 
case) or ~ 0.7uJy (radio-quiet case), ^ 22 mag (AB), 
and 6.4 x 10"^^ ~ 3.6 x lO^^^ ergs cm^^ for 1.4GHz, 
24/xm, and 2-8keV, respectively. The Icr detection limits 
for the available data at these wavelengths, however, are 
8.5uJy, 21.2 mag (AB), and 6.7x lO'^^ ergs cm-^ g-^ for 
1.4GIIz, 24/xm, and 2-8keV, respectively. Non-detections 
of TENIS-ZDl at these wavelengths therefore cannot rule 
out a dusty radio-quiet type-2 AGN with weak X-ray lu- 
minosity at z ^ 1.5. Hence we did a further check by 
doing an SED fitting within a redshift range between 1.0 
and 2.0 for TENIS-ZDl us ing an AGN/quasar template 
from iPolletta et al. I (|2007f ). The best-fitted template is 
the QS02 template. According to the fitting result, the 
expected z-band magnitude is ^ 26 mag, which should 
be a 5a detection in the GEMS F850LP image, but the 
signal-to-noise ratio of TENIS-ZDl in the GEMS im- 
age is less than la. The best-fitted SED also fails to 
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Fig. 10.— The best-fitted SED and the x vs. redshift plot. The 
black data points indicate the AB magnitudes of TENIS-ZDl in 
16 bands. The red, green, and blue lines indicate the best-fitted 
templates of galaxy, dwarf star, and AGN/quasar, respectively. 
The plot for the best-fitted galaxy template shows that the 
photometric redshift fitting of TENIS-ZDl has a very low from 
^phot = 7.0 to Zphot = 9.0, while the weighted ^p^ot 7.8. 

photometric errors. We therefore discuss the possibil- 
ity of that TENIS-ZDl is a low-z galaxy. The red data 
point in Figure [S] shows the position of TENIS-ZDl in 
the J — 3.6/im vs. F850LP — J diagram. According to 
Figure [H although TENIS-ZDl is right on the bound- 
ary of our selection criteria, it is actually a lower-limit 
of its F850LP — J color since it has no detection in the 
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reproduce the steep color slope b etween Kg band and 
4.5/im. Furthermore, Figure 13 in ICaoak et all ([20lTh 
shows that objects at z ^ 7 still can match the z > 7 
selection criteria. The J — 4.5/im color of TENIS-ZDl is 
1.33 ±0.24, which is outside the color range of any z ^ 7 
galaxies (J — 4.5/im < 1.0 even for extreme cases). It 
is therefore extremely difficult to explain TENIS-ZDl 's 
photometric properties using low-z galaxies. 

4) Galactic cool stars: Although our selection criteria 
can separate dwarf stars and high-z galaxies (see Fig- 
ure [6]), contaminations from Galactic cool stars can still 
be present because of photometric errors. The red data 
point in Figure [5] indicates TENIS-ZDl. According to 
Figure H the 3.6/im - 5.8/tm color of TENIS-ZDl is 
too red as compared to that of dwarfs. Moreover, the 
FWHM of TENIS-ZDl in the J-band image shows it 
is an extended source. All these evidences suggest that 
TENIS-ZDl is not a Galactic cool star. 

5. DISCUSSIONS 

5.1. Surface density 

The cumulative surface density of our sample is ^ 
1 X IQ" ^ per arcm i n^ to J < 25.1. According to Fig- 
ure 4 in lYan et al. I ()2011|) (Yanll, hereafter), this value 
matches to that estimated from the LF based on the 
WFC3 Z85o-dropout results for z ~ 7 galaxies. Yanll 
claims that the estimation of the cumulative surface den- 
sity is applicable over the redshift range of 6.4 < z < 7.7. 
In other words, it suggests that such an extremely lumi- 
nous galaxy with 6.4 < z < 7.7 can be found in a ^ 1000 
arcmin^ survey. Therefore, the discovery of TENIS-ZDl 
is statistically predictable based on previous studies us- 
ing fainter samples. 

5.2. UV luminosity and star formation rate 

EAZY code cannot provide physical properties (e.g., 
age, stellar mass, etc.) other than photometric redshifts. 
Therefore we use the New- Hyperz (Rose et al. B) w ith 
the GALAXEV templates (|Bruzual fc Chariot l[200l to 
estimate the physical properties of TENIS-ZDl. We pass 
the best-fitted ZpUot derived by EAZY to New-Hyperz 
and treat it as a fixed parameter during the template fit - 
ting. The extinction model from Calze ttTet al. I ()2000l ) 
is used. The best-fitted template {y^ = 1.38) shows 
that TENIS-ZDl is a starburst galaxy, with an age of 
only 45M years and a stellar mass of 3.2 x 10^° Mq. 
The extinction is mild, with an Av = 0.6. The absolute 
UV magnitude Migoo is -22.35 before the extinction cor- 
rection. We then derive the extinction-corrected Mieoo, 
MJgQQ, using the following formula: 

Fobs =FoX l0i-OAxAvxK,/4.O5) (3) 

(jCalzetti et al. Il2000( ). where Kx = i^ieoo = 9.97, and 
Av = 0.6. The correction value is -1.48 mag and hence 
M?[g{i5 is -23.83. We then estimate the SFR using the fol- 
lowing formula (jMadau. Pozzetti. fc Dickinson Ill998| ): 

SFR = Luvierg s'^ Hz-'^)/{8 x lO^^^). (4) 

For the extinction-corrected case, the estimated SFR is 
^ 200 M0 year~^, while for the extinction-uncorrected 
case, the SFR is ~ 50 Mq year"^. 

^ |http : //www ■ ast ■ obs-mip ■ f r/users/roser /hyperz/ 1 



We try to estimate the number of galaxies at z>7 us- 
ing th e semi-analytic simulation result from iGuo et al. I 
(|2011f) . which is Ijased on the Millennium Simulations. 
The result suggests that ONE galaxy at 7<z<8 with 
a stellar mass greater than 10^" M© can be detected 
in the TENIS project, which supports the existence of 
TENIS-ZDl. However, the value of as used in the Mil- 
lennium Simulations is 0.9, and the most recent estimate 
of gs is 0.809 ± 0.024 derived using the WMAP7 data 
(jKomatsu et al. 11201 ll ). which is almost 4cr smaller than 
that used in the Millennium Simulations. Structure for- 
mation would be faster with a larger cg as massive ob- 
jects can form earlier. It is possible that, with the revised 
(jg, the number of galaxies similar to TENIS-ZDl pro- 
duced in the Millennium Simulation might be less than 
one in our survey volume. 

There are several massi ve galaxies at z ~ 6 with 

IRAC detections re ported in lYan et al.1 (|2005l [20061 ) and 
lEvles et al. I (|2007f ) . According to their analyses, the stel- 
lar masses of these galaxies are about several 10^° Mq 
and their progenitors can be observed at z > 7. In par- 
ticular, they point out that the z ~ 6 sample strongly 
indicates that the universe was already forming galaxies 
as massive as ~ 10^° Mq at z > 7 and possibly even at 
z ~ 20, which supports the existence of TENIS-ZDl. 

Although it is challenging to explain how to accumu- 
late so much materials and trigger such violent star for- 
mation at z > 7 with current cosmological models, the 
observational results of z ~ 6 galaxies suggest the ex- 
istence of TENIS-ZDl is reasonable. Therefore bright 
high-z objects like TENIS-ZDl can provide an impor- 
tant constraint for developing and modifying models. 

5.3. Luminosity function 

Based on the value of Migoo, TENIS-ZDl is an ex- 
tremely luminous object at z > 7. To date, TENIS-ZDl 
is brighter than all p ublished z > 7 samples except for the 
sample provided bv lCapak et al. I (|2011[ )fCapakll. here- 
after). Hence TENIS-ZDl provides a new constraint to 
the very bright end of the LF at z > 7. In order to make 
an LF plot with TENIS-ZDl, we estimate the effective 
comoving volume for our analyses using a Monte-Carlo 
simulation which is modified from Ouchi09. 

First of all, we generate a mock catalogue for Lyman- 
Break Galaxies (LBGs) usi ng the starburst templ ate of 
the GALAXEV model f Br uzual fc Chariot I [2001 with 
apparent J magnitudes from 22.0 to 27.0. The surface 
density vs. apparent magnitude relation of the mock cat- 
alogue follows that of all the objects shown in Figure [2l 
where we have assumed that LBGs have a similar sur- 
face density vs. apparent magnitude dependency. The 
objects are uniformly distributed in the redshift space, 
and their ages are randomly distr ibuted from IM years 
to 500M years. The extinctions of iCalzetti et al. 1 pOOOl ) 
with Av = 0.0 to 2. are also applied to s i mulate bluer 
and redder LBGs. iDow-Hveelund et al.1 (|2007f) show 
that roughly only one-third of dropout-selected galax- 
ies at z = 5 — 6 have strong Lya lines with rest-frame 
equivalent widths > 20A. The fraction might be even 
lower at z > 7 because of the absorption by IGM. We 
therefore added Lya emissions with equivalent widths of 
20A for 25% of the objects. We assume that 50% of the 
objects have no Lya emission. For the remaining 25% 
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Fig. 11. — Luminosity function at z ~ 7. The red filled 
square indicates our sample, the black filled square indicates 
the sample in Ouchi09, and the black cross indicates the sam- 
ple in Capakll. The other gray data points sh ow the sam - 
ples from Bouwcns ct al. (2008, 2010b, c); McLure cfaQ ll2010l') : 
lOesch ct al. (2009, 2010a, b); CastcUano ct al. (201(|). The blue 
solid line and green solid line represent the LF at 2 ~ 7 pro- 
vided by Yanll and OuchiOQ respectively. The LF at z ~ 6 in 
IBouwens et al. I II2007I ) is also shown as the gray solid line. Our 
result is consisitent with both OuchiOQ and Yanll. 

of the objects, we added Lya emissions with equivalent 
widths of 1-20A. After taking aU of the abovementioned 
effects into account, we disturb the photometry of the 
catalogue to match to the real data. 

Equations[T]and[2]are applied to the mock catalogue to 
select candidates. We computed the redshift distribution 
C{m,z) as the ratio between the number of candidates 
and the input objects in a magnitude bin of m = 251q25, 
and estimated the comoving volume density for a given 
magnitude bin by using the relation 



n{m) 



Ncandjm) X Fc{m) 



(5) 



which is similar to Equation 2 in Ouchi09, where 
Ncand{ni) is the number of candidates in the given mag- 
nitude bin, Fc{m) is the completeness correction, ^ is 
the differential comoving volume for the field size of the 
ECDFS, and m = 25 is the only magnitude bin we used 
since our candidate is in that magnitude bin. In order 
to convert the unit of n(m) {N 0.5mag~^ Mpc^^) to 
N mag~^ Mpc^^, a factor of 2 needs to be applied to 
n{m). We make the LF plot, including the data from 
several previous studies, as shown in Figure [TT] The LFs 
provided by Ouchi09 and Yanll fitted including lumi- 
nous samples are also plotted. The error of our data 
point is calculated based on a simple Poissonian uncer- 
tainty for an event rate of 1, which is ±1. On the other 
hand, the 90% Baysian c onfidence level for such an event 
rate (jKraft. Burrows, fc Nousek 1991) would be 0.08- 
3.93. In either case, our result is consistent with both 
Ouchi09 and Yanll. 

Several papers claim that the LF at z ^ 
7 de creases from z ~ 6 at the brigh t end 
(e.g., iMannucci et al. I 120071: iCastellano et al. I 120101: 



lOuchi et al. 112001 IBouwens et ani2010bD . while Yanll 
shows no/ very mild evolution. These studies are done 
using samples with Mjjv > ~22 while our sample 
would provide a better constraint at the very bright-end 
{Mjjv < —22). However, our sample is not sufficient 
to distinguish between the two cases since our result is 
consistent with both Ouchi09 and Yanll within the un- 
certainties. On the other hand, Capakll shows a z ^ 7 
LF with a very high bright end, according to a sam- 
ple even brighter than ours. We note that we use very 
restricted criteria to select our candidate. While these 
criteria would eliminate contaminations as much as pos- 
sible, some objects at z > 7 could also be rejected at the 
same time. From this point of view, the constraint that 
our sample provides for the LF is a lower limit, and thus 
we cannot completely rule out the LF from Capakll. 
Furthermore, cosmic variance can be as high as a factor 
of 10 since the ext remely luminous galaxies are strongly 
clust ered at z > 5 (jOuchi et al. 1200 j : lHildebrandt et all 
|2009( ). It is reasonable that the results based on observa- 
tions covering less than a few square degrees may show 
large discrepancies. 

6. CONCLUSIONS 

The TENIS project provides deep J and Kg images 
for the ECDFS. We use these data with the other deep 
surveys in optical and infrared to find z > 7 objects. 
New color criteria with IRAC data for selecting z > 7 
candidates are used because we find that accurate IRAC 
photometry provide the key to resolve contamination by 
\ow-z galaxies and Galactic dwarf stars. Because the 
PSFs of IRAC data are much larger than those of the 
data in the other wavelengths, we introduce a novel de- 
convolution method to address the confusion in the IRAC 
images and provide accurate IRAC fluxes. After carefully 
checking our sample, we found one candidate at z > 7, 
TENIS-ZDl. 

The weighted Zphot of TENIS-ZDl is 7.8, with an esti- 
mated stellar mass = 3.2 x 10^° M©. The extinction- 
corrected current SFR of TENIS-ZDl is ~ 200 M© 
year~^, but must have been very much higher in the 
past. We summarize the results of our sample below: 1) 
The discovery of TENIS-ZDl is predictable in an 1,000 
arcmin^ survey according to the surface density derived 
from literatures. 2) Our sample matches to both LFs 
from Ouchi09 and Yanll. 3) While the existence of 
TENIS-ZDl is supported by the observations at z ~ 6, 
it is still hard to explain how such a massive galaxy can 
form in the early universe based on current cosmological 
models for structure formation. 

In the end, we would like to emphasize that the TENIS 
project provides the most comprehensive dataset possible 
to address the problem of flnding luminous high-z galax- 
ies. Compared with fainter high-z samples, such high- 
z massive objects like TENIS-ZDl provide the greatest 
leverage in testing galaxy formation/evolution models. 
Such objects are the best candidates for follow-up spec- 
troscopy with large telescopes; at the same time, to min- 
imize wastage of observing times on large telescopes, all 
efforts must first be made to weed out contaminants. The 
sizes of the z ~ 7 bright samples are very small due to 
the limited field size of deep near-infrared surveys, and 
the purity and completeness issues of the dropout tech- 
nique need to be further investigated. Due to the severe 
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contamination issue of the bright end of the 2: ~ 7 LF, 
all the studies make great efforts on increasing the pu- 
rity of their samples, however, all these efforts may also 
harm the completeness of the samples, which makes the 
bright end of the z ~ 7 LF seriously under-estimated. 
iHu h Cowie I ()2006l ) point out that the bright end of 
the LF at 2; ~ 6 based on LAE studies is much higher 
than that based on LBG studies. Capakll also claims 
that the bright end of the LF at z ~ 7 could be much 
higher than the results derived using the dropout tech- 
nique. The bright samples will increase very soon since 
there will be more and more large field near-infrared sur- 
veys. However, simply increasing the bright sample size 
cannot solve the bias issue. The best way is to confirm 
their redshifts spectroscopically. Each spectroscopically 
confirmed object will set a robust lower-limit for that lu- 
minosity bin of the LF. The James Webb Space Telescope 
(JWST) and 30m-class telescopes (e.g., TMT, GMT) will 
be certainly capable for these confirmation observations. 
But before the era of JWST and 30m-class telescopes 
comes, the current lOm-class telescopes are still able to 
perform the observations for objects with Mjjy < —21 
as long as strong Lya lines exist. These spectroscopic 
observations will provide very important information for 
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